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Abstract 

We develop the Redfield equation for delta-correlated gaussian noise and 
apply it to the case of two neutrino flavor or spin precession in the presence of 
a noisy matter density or magnetic field, respectively. The criteria under which 
physical fluctuations can be well approximated by the delta-correlated gaussian 
noise for the above cases are examined. Current limits on the possible neutrino 
magnetic moment and solar magnetic field suggest that a reasonably noisy 
solar magnetic field would not appreciably affect the solar electron neutrino 
fiux. However, if the solar electron density has fluctuations of a few percent 
of the local density and a small enough correlation length, the MSW effect is 
suppressed for a range of parameters. 

1 INTRODUCTION 



Neutrino oscillations in the presence of matter and magnetic fields have been an 
area of intense study for approximately the last ten years. In the Mikheyev - Smirnov 
- Wolfenstein (MSW) effect, electron neutrinos on their journey from the core, are 
resonantly transformed into muon or tau neutrinos H]. If neutrinos are Majorana 
fermions with transition magnetic moments, they can undergo a magnetic resonant 
transformation into muon or tau antineutrinos 0]. Neutrinos from a supernovae 
explosion can be transformed from one flavor to another as they pass through the 
outer part of the star [Q. In many stellar situations, the matter density and/or 
magnetic fields may fluctuate about a mean value. Some well known examples where 
fluctuations are likely to exist include the magnetic field and the matter density in the 
solar convective zone and also the turbulence of the post-supernovae matter which 
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has been blown off by the explosion. A general approach to the neutrino oscillations 
in inhomogeneous matter was developed in Ref. 0. A Study of matter fluctuations 
which are not random, but harmonic 0], or occur as a jump-like change in the solar 
density [||, are available in the literature. Matter currents and density changes effect 
neutrino flavor oscillations in a similar way and have also been examined in Ref. 
Although matter current effects become important only if the velocity is somewhat 
close to the speed of light, noisy mixing of matter would also mimic a fluctuating 
matter density. A priori, fluctuations in such fields may be well approximated by 
random noise added to an average value. In this paper we show how such noise will 
affect neutrino oscillations for the situation in which the correlation length of the 
randomly fluctuating part of either the matter density or the magnetic field is small 
compared with the neutrino oscillation length. 

The case of neutrino spin precession in a noisy magnetic field was considered by 
Nicolaidis for neutrinos in vacuum |p. The noise was taken to be well approximated 
by a delta-correlated gaussian distribution with the result that the normal oscillations 
become damped with a relaxation time of tj-ei = < > Tc)~^, where Br is the 
randomly fluctuating part of the magnetic field and Tc is its correlation time. Enqvist 
and Semikoz considered neutrino oscillations in randomly fluctuating magnetic field, 
approximately, by averaging the coefficients in the third order differential equation 
for the z-component of the neutrino spin, in a constant matter density [l^. Their 
result suggested that the effect of the randomly fluctuating part of the magnetic field 
was similar to the effect of a larger (constant) matter density, hence reducing the 
neutrino precession. Our result is quite different; namely that both a noisy magnetic 
field and a noisy density (for a constant averaged density) act to depolarize the 
neutrinos. For a noisy magnetic field, if the probability of transition is greater than 
one half without the random fluctuations, the inclusion of the random fluctuations 
will reduce the transition probability, and if the transition probability without the 
random fluctuations is less than one half, the inclusion of the fluctuations will increase 
the transition probability. If the randomly fluctuating part of the magnetic field is 
strong enough or is allowed to act for a long enough time, a complete depolarization of 
particles-antiparticles will occur. For a noisy density, we find that the MSW transition 
probability is suppressed. For the case of strongly adiabatic MSW transitions and 
large fluctuations, the averaged transition probability saturates at one half. 

In Section II, we develop the equations governing the time evolution of the av- 
eraged probabilities of being found in the Nth level of an N-level system, subject to 
random and non-random potentials. These equations are derived for the case where 
the randomly fluctuating part of the field is taken to be a delta-correlated gaussian. 
In Section III, we show several analytically solvable examples for the case of a two 
level problem, and numerically examine cases the analytical solutions of which are 
not instructive. In Section IV, we investigate the conditions under which a real, 
physical fluctuating field will be well approximated by the equations developed for 
the delta-correlated gaussian case. We then apply these conditions and examine the 



cases of spin-flavor and flavor precession of neutrinos in the sun. Section V presents 
a discussion of results and our conclusions. 



2 FORMULATION OF THE APPROACH 



We first consider the general case of an N-level system with time dependent level 
splitting and transition terms, in the presence of an additional randomly fiuctuating 
term. Since we are interested in the ensemble average of the probability of finding 
the system in a given level at time t, we consider the density operator defined as, 
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where. 
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Assuming that ip obeys the Schrodinger-like equation. 



p obeys the equation. 



where the Hamiltonian is taken to be linear in fluctuating fleld, i.e. 

H ^ Hoit) + B{t)M' , 
where M' is independent of time. Eq. (4) becomes. 
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where. 



and C/q satisfies. 



i-^Pi = [Hi, Pi], 



Pi = Ul p Uo, 
M = f/(J M' Uq, 
Hi = B{t)M, 

ijlJo = Hollo. 
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dhdt2B{U)B{t2) [M(ti), [M(t2), Po]] + ■ • • (9) 



This equation can be solved by iteration as, 

pit) = po-t J' dhB{t^)[M{t^),po] 

- IT 
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We now assume that B(t) is such that, 

<B{t)>=0, <B{h)B{t2)>=a^fi\t2~h\), (10) 

with the average of all higher odd products of B vanishing, and all higher even prod- 
ucts given by the sum of all possible products of f's for example, 

< B{h)B{t2)B{h)B{U) >= a^lfufsA + /13/24 + /14/32], (11) 
where, /12 = /(I ^2 — I)- This will result in only the even products contributing to 

< m >. 

Further analytic progress can be made if the random potential is a delta correlated 
gaussian distribution, where /(| ti — t2 |) becomes 2r(5(ti — t2). Averaged values of 
functions of the random field B(t) can then be expressed as a path integral, 

V[B{tMB)e-^o'^^^ (12) 

-00 

where, 

V[Bit)]=Y[dBiU)J^, (13) 

i ' 

and k = IcP't. This reproduces the form of the averages of the even power products, 
and the integrals in Eq. (9) can be explicitly calculated with the result, 

<pi{T)>=po - a^T [ dti[Miti),[Mih),po]] 





+ rdt,[M{h), [M(ti), r dt2[M{t2), [M{t2),po]]]] (14) 
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This, however, is just the iterative expansion of the differential equation, 

^ < pj{t) >= -aV[M(t), [M(t), < pi{t) >]]. (15) 
One may equally well express this result as, 

^ < m >= ~a'T[M\ [M', < m >]] - t[Ho{t), < m >]. (le) 

Eq. (16), in general, defines a set of coupled first order linear differential equations 
with the constraint that Trp = 1. This is a matrix form of the Redfield equation |jlT| 



In the case of neutrinos moving in a varying background, Eq. (16) is equivalent to 
Eq. (49) in Ref. §. 



3 Analytic Examples 



Perhaps the simplest examples are those for which one could obtain an integral 
expression for the exact solution for the probabilities of the system being found in 
the N*^ level at time t, and then perform the averaging using Eq. (12). Such an 
expression is possible, in general, only in the case where all N levels are degenerate. 
For example, consider the case in which the Hamiltonian is an N x N matrix which 
is of the form H{t) = {Bo{t) + Br{t)) x M, where M has zeros on the diagonal 
and ones everywhere else. One could, for this case, obtain an integral expression for 
the probabilities and then average them. Instead, we take i^o(^) = -Bo(t)M giving 
^o(^) = 6xp(/q (iti?o(t)M) and use Eq. (15) which becomes, 

I < Pi{t) >= -'^[M, [M, < m >]], (17) 

where 

k{t)/2 = a^{t)T. (18) 

Taking the commutator of M with ^ < pi{t) >, and noting that — [N — 1) -\- 
(N - 2)M, one obtains, 

|[M, < pj{t) >] = -^k{t)[M, < p:{t) >]. (19) 
This equation is easily solved with the solution, 

< p{T) >= X{T) - Ar-2(i _ e-^/o''^*'=(*))[M, [M,X{T)\\ (20) 

where X{t) — Uo{t)p{0)UQ{t). This reproduces the result of Nicolaidis for the case 
where M' is the Pauh matrix ax, and the initial state is ■01 (0) = 1, '02(0) = 0. 

A second analytically solvable example is a two level system, with an arbitrary 
time dependent level splitting, in a purely random field. This would correspond to 
the case of matter-enhanced neutrino spin precession in a noisy magnetic field. In 
this case, there exists no analytic expression for the probabilities which one could 
in principle average. We begin by using Eq. (16) with Hq — A{t)az + Bo{t)ax 
and M' = a^. Defining r{t) = |(< pn > — < P22 >),x — 2Re{< pu >) and 
y = 2Im{< pi2 >), one obtains. 
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If So = Eq. (21) is trivially solved giving. 



r(r) =r(0)e-'/o^'^*'=(*), 



(22) 



resulting in, 

< \MT)\^ >= -(1 - e-'i'o''^*'^^*)) + I V'l (0) I V'io^ (23) 

This result is surprising in that it is independent of the level splitting and identical 
to the result for zero level sphtting. If Bq ^ 0, one could still obtain an analytic 
expression for the averaged probabilities if the level splitting is constant, although 
such an expression is not instructive due to the complicated nature of third roots. 
We instead calculate numerically cases of constant level splitting. Figure (la) through 
(Id) show the probability |-?/'i(t)p for A = 0, A = 0.5, A = 1 and A = 10 respectively, 
for Bq — \ and k — 0.05. One observes from these figures that as A is increased, 
the oscillations become heavily damped while the exponentially decreasing upper 
envelope remains essentially unchanged. This qualitatively indicates that even when 
one is far from resonance (A = for the case in which A is time dependent), one can 
still obtain complete depolarization for such a randomly fluctuating fleld. 

Another possibility occurs if one has the case of a noisy diagonal term with a non- 
noisy off-diagonal term, namely H = (Aoit) + Ar{t))az + BQ{t)ax- This case would 
correspond to a noisy matter density in the context of matter-enhanced neutrino 
oscillations. Eq. (16) becomes. 
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If Ao{t) = and Bq is time independent, one can obtain an exact solution which, for 
the initial condition t/'i(0) = 1,1^2(0) = and 4Bq > k"^ is, 

11 k 

|^i(i)|2 = _ + -e-^\-s\nujt + cos ujt], (25) 
2 2 oj 



where LO = y4:BQ — k'^. We again numerically calculate examples in which the non- 
random diagonal term is non-zero. Figures (2a) through (2d) show |-?/'i(t)p for Aq — 
0,74o = 0.5, = 2 and Aq — 10 respectively for Bq — 1 and k — 0.05. In in 
contrast to the case of a noisy off-diagonal term, in this case one observes that as 
the oscillations are damped, the exponential relaxation time increases. Thus, if one 
is far from resonance, the depolarization is very much suppressed, as is physically 
reasonable. 



4 NEUTRINOS IN NOISY SOLAR FIELDS 

In Appendix I, we detail the expansion of the density matrix in the case where 
the randomly fluctuating fleld has a flnite correlation time, with the result that Eq. 
(16) can accurately describe a real physical situation if, 

{d/dt{\oga{t)) + [^o(i), M'])Te < 1, (26) 



for all t of interest, where is the correlation time of the randomly fluctuating field 
and a is it's root mean square value, which may be time dependent. If this is satisfied 
for a given problem, is given by, 

k{t)=(5o?T,{t), (27) 

where /? is a factor of order unity, equal to 2 for our choice in Appendix I of a theta 
function correlation. 



4.1 Spin Precession in a noisy solar magnetic field 

We first consider the case of matter-enhanced spin-flavor precession of Majorana 
neutrinos with a negligibly small MSW mixing angle. In this case the evolution 
equation is two by two, and is given by H, 




Am 
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where Am^ = m2 — mf with rrii the masses of the mass eigenstate neutrinos, E is 
the neutrino energy, /i is the neutrino transition magnetic moment, B and Br are the 
magnetic field and its noise, and Qe and are the matter potentials given by, 

ae = -^Gi.(2iVe-iV„), a^ = zlG'^iv„, (29) 

where Nf.{r) and Nn{i) are the electron and neutron number densities. 

From the discussion in the previous section of a noisy off-diagonal term, we can es- 
timate whether a reasonably noisy magnetic field can effect the transition probability. 
In order for the random fluctuations to have much effect, 2k Ar ~ 1. The condition 
for this real process to be approximated by a delta-correlated gaussian, namely Eq. 
(26), results, for a constant rms value of the the magnetic field, in the condition, 

A 2 

Tc X [-^ + V2Gf{N, - iV„)] « 1. (30) 
Perhaps the most likely place a noisy magnetic field would exist is in the solar 



convective zone which extends from approximately 0.7 Rq to the surface [|T2|. The 
magnitude of the solar magnetic field observed at the surface of the sun can reach 
local values of several kilogauss in magnetic storms, and it is thought that it may 



reach values of 100 kG near the bottom of the convective zone |T^. While a large 



lOOkG ) field which extends through out a large fraction of the convective zone 



is not ruled out, such a large, extensive field is thought to be unlikely [|T4|. Field 
strengths as large as several times lOOkG could exist if limited to extensions of about 
10 — 100 km. In order to estimate the maximum effect of magnetic field fluctuations. 



we will assume that the solar magnetic field from 0.7Rq to O.85i?0 can be considered 
to be randomly fiuctuating with zero mean and a rms value of lOOkG. We take the 
neutrino transition magnetic moment to be = 3 x 10^ ^^yU^, in accord with the 
maximum bound from plasmon decay in pre-helium fiash red giant stars ||15|. By 
choosing a value of A'm'^/2E, one can determine the maximum Tc (assumed to be 
constant) such that Eq. (30) is satisfied throughout the above region in the solar 
convective zone. Therefore, we take Tc to given by, 

r, = 0.1 X (I - ^ + V2GFiN, - Nn)\ma.r'- (31) 

Given that Eq. (30) is obeyed, the average electron neutrino probability after travers- 
ing this region is given by Eq. (23), with ipei^) = 1, and for small 2k Ar the muon 
antineutrino average probability is approximately given by kAr. In Figure (3) we 
plot the value of kAr = 4 < {^BY > Tc{0.15Rq) as a function of Am? /E. One 
observes that the maximum occurs at about 1 x 10~^eV^/MeV corresponding to the 
largest permissible value of Tc. The peak is due to the fact that the permissible values 
of Tc are maximized when the neutrino is in resonance somewhere between 0.7 Rq and 
0.85Rq. If one wished to satisfy the condition on Tc for Am? / E less than or equal to 
4 X 10~^eV^/MeV, one would obtain about 3 % muon antineutrinos on average. This 
is the same order of magnitude as that is produced by resonant transitions in this 



region of the sun, for comparable magnetic moments and magnetic fields pTS] , 0, |l4|. 

We would also like to point out that in the case of non-zero mixing between 
neutrino species in a purely random magnetic field like that considered above, one 
can derive the exact (assuming Eq. (26) holds) time dependence for the sum of the 
neutrino probabilities i.e. |'?/'e(r)p + \il)^{r)\'^. The Hamiltonian is given by [^, 

H-i",^ l)^.B,i L f V (32) 



where, M' = icTy and. 



If one evaluates Eq. (16) one can derive the equation for the sum of the neutrino 
probabilities. The solution is given by Eq. (23), with < |?/'i(T)p > replaced by the 
sum of the electron and muon neutrino probabilities. Therefore, the same estimates 
as in Fig. (3) apply to this case but for the sum of the antineutrino averaged prob- 
abilities. One should note that the fiavor precession will, in general, be changed by 
the inclusion of the magnetic field. 



4.2 Flavor Precession with a Noisy Matter Density 



The equation governing matter-enhanced flavor oscillations (for two flavors) is 
given by, 

-^<^os29 + ^^G,K ^sm20 \ ( ^.A .... 

dr U.J 1^ ^sin2^ ^cos2^ - ^G^TVe j J ' 

where 6 is the vacuum mixing angle. We investigate the case in which A^e = (1 + 
P)Ne{r) where /? is a random quantity whose rms value is the rms value of the density 
fluctuations relative to the average density and A^e is the electron density in the 
standard solar model. We satisfy Eq. (26) by choosing, 

n = Q.lx{^sm29)-\ (35) 

where the logarithmic derivative in Eq. (26) is much smaller for the sun than the 
term we have kept. We numerically solve Eq. (24) for the electron neutrino survival 
probability (averaged for infinite distance beyond the surface of the sun), as a function 
of Am^/ E for fixed values of sin^ 29 and < j3'^ >, with r^. given by Eq. (35). We have 
used the solar electron density of Bahcall and collaborators. Figures (4a) through 
(4d) show the electron neutrino survival probability as a function of Arn^/E for 
sm^29 = 0.8,0.1,0.01 and 0.001 respectively, with V< > = 0.02 In these figures, 
the curve which reaches the smallest values in the center of the plot (from about 
Am'^/E = 1 X 10"^ to 1 X 10""'^cV^/MeV) is the probability in the absence of the 
density fiuctuations. One observes a suppression of the fiavor transition in this region 
by as much as 20 % in Fig. (4c). The region of the largest effect (around Aw? /E — 
1 X 10~^ is in agreement with the maximum value of 2kAr, where Ar is the width of 
the resonance region. In order to show the sensitivity of the suppression to the size 
of the fiuctuations, we show the same plots as Fig. (4a)-(4d) in Fig. (5a)-(5d), with 
the exception that the rms fiuctuation in the density is doubled to 4 % of the average 
density. One observes a very large effect in the previously mentioned region, which 
saturates at a probability of one half, for values near Avn? jE = 1 x 10~^eV^/MeV. 
One should note that in Figures (4a)-(5d), Tc was chosen for each value of Avn? jE to 
obey Eq. (35) . Thus, the chosen for Ain^ / E = Ix 10~^eV^/MeV, while permissible 
for smaller values of Am'^/E, would not have given as large an effect, and would, for 
larger values of Aw? /E violate Eq. (35). 

In Figs. (4) and (5), there appears to be an enhancement of the transition prob- 
ability for values of Aw? jE > 1.5 x 10~^. Neutrinos of with these parameter values 
do not go through a resonant transition in the sun, since a resonant transition would 
require a density larger than the maximum solar density. However, to a varying ex- 
tent they are within the second half of their resonant region near the center of the 
sun. The solar density profile near the center is rather flat and therefore one obtains 
an enhancement of the transition probability similar to that in Fig. (3c) and (3d). 



One needs to question whether such fluctuations can arise in the sun, and whether 
the correlation lengths, which can be quite small (r^ for Am? /E = 1 x 10~^eV^/MeV 
and sin^26' = 0.01 is about 10km ) could be realistic. If one were to do a numerical 
study utilizing a good random number generator, one would not be limited by Eq. 
(35) and could consider any correlation length. Our purpose is only to show that such 
effects could he important. 



5 Conclusions 

We have derived a Redfield differential equation for the time dependence of av- 
eraged values of functions of the probability of finding an N level system in the N*'^ 
level after being subjected to a randomly fluctuating field, of the delta-correlated 
gaussian type. This formalism applies as an approximation to the case where the 
fluctuating field can be described by a finite correlation time. This approximation is 
valid if the product of the correlation time and the energy scale of the Hamiltonian is 
small. Upon applying this to neutrino flavor or spin-flavor precession, we have shown 
that the probability will relax eventually to a value of one half, if the neutrino spends 
a sufficiently long time in a medium with randomly fluctuating matter density or a 
randomly fluctuating magnetic field. In the case of a fluctuating magnetic field, the 
relaxation is independent of the matter density, assuming the correlation time of the 
field fluctuations is small compared to the neutrino oscillation length in matter. In 
the case of fluctuations added to a constant matter density, the probability again will 
eventually relax to a value of one half, but the relaxation time is greatly increased if 
one is far from the resonant condition. 

We have also examined the case of neutrino spin-flavor precession in the sun, for 
a purely random magnetic field and no flavor mixing. It appears that the current 
limits on the neutrino magnetic moment and the guesses concerning the maximum 
values of the solar magnetic field, combine to give only a small effect on the aver- 
age electron neutrino flux. When there is flavor mixing and a purely random field, 
the combined average probabilities of the neutrinos is again a simple exponentially 
decreasing function and therefore the results of a purely random field and no flavor 
mixing apply to this case as well. For the case of a randomly fluctuating electron 
density, the MSW effect can be strongly suppressed for rms fluctuations of 4 % of the 
local electron density. However, this requires correlation lengths of about 40km, and 
seems to give a significant effect only for neutrinos which have their MSW resonant 
transition deep in the sun. 

In spite of these problems, we believe that a numerical study, which should give 
similar results for correlation times which do not badly violate Eq. (35), may bring 
to hght many interesting effects. Implications of the density fluctuations discussed 



here on stellar collapse and supernova dynamics will be published elsewhere |T7|. 
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Appendix 



For -^o(^) given in Eq. (5), one selects the element of largest value during the time 
the system is in the presence of the random field. Let this largest value be Ejjiax- We 
consider the case where rEmax <S 1. We take fij in Eqs. (10) and (11) to be given 

by, 

hj = OiT-\U-tj\), (36) 

where again the average of odd products vanish. We rewrite the average of Eq. (9) 
as, 

< Piit) >= po+ < pf\t) > + < pf\t) > + <pf\t)>+--- (37) 

where, 

Jo Jo 

+ r dtiF2N[M{t2N): [M{t2N-l): ' ' ' M (t,)] ■ ■ ■]] (38) 

Jo 

and, 

P(1,2,---,27V) 

F2N — fmnQfnsnifnsne ' ' ' fn^N-iniNi i"^^) 

ni---n2N 

where, P(l,2,- • • ,2N) means all permutations. We explicitly show the second and 
third terms: 

< pf\t) > - -a^ tdh[M{h), r [M{h) + ^^^{t2 - h),po]] 
Jo Jti-T at 

~ -aV l^dh[M{h),[M(T,),po]] + a'^ dh[M{h),[^^^, po]], (40) 



where dM(ti)/dt = iU^[Ho, M']f/o and is therefore proportional to Emax- In the third 
term one has a sum of three products of two f 's, (/12/34 + /13/24 + /14/32) only the 
first of which has the times in the order of the times appearing in the nested integrals. 



The first of tliese tliree terms gives, 



<P^'\t) >12,34 ~ aV fdh[M{h),{[M{h), r dt2[M{t2),[M{t2),Po]]]] 

JO Jo 
~ l(l^^J^'dHmt2),[M{t2),Po]]]] 

+ [M{t^),l^'dUM{h),[^^,Po]]]] (41) 
+[Mih), [M{h), [Mih),po]]]])} + O(r^). 

Tlie last term, liowever, is not proportional to Emax- This will cause a problem since 
it could contribute to the next order term of the average . It turns out that it cancels 
the largest term coming from the remaining two f products. That these remaining f- 
products arc of largest order can be seen by noting that when the argument of the 
theta function in an f connects two non-sequential times, the intermediate time(s) 
must also be within r of the larger time in the theta function. For example, the 
remaining two contributions to < pf\t) > are, 

< pf\t) >13,24 + < pf\t) >14,23= 

Jo 

X / dt2 dtJ dt4+ dt4) (42) 

Jtl-T Jti-T ^Jh-T Jti-r ' 

= Oi^% fdh[M{h),[M{t,),[M{h),[M{h),po]]]]. 
2 Jo 

This feature appears to continue throughout each term in the entire expression. 
Therefore, if rEmax <S 1, and neglecting terms of order rEmax and smaller, 

<pf''\t)> = a'"" t"" f dhNm^N) AM {hN). 

Jo 

rt2N ^ ^ ri2N-2 

/ dt2N-2[M{t2N-2,[M{t2N-2), ■■■ (43) 

Jo Jo 

rdt2[M{t2)AM{t2) Jo]] ■■■]]]]■ 

Jo 

which leads to Eqs. (15) & (16). 

In the above derivation we have assumed a to be a constant. If a depends on 
time, the condition, tE^^x <S 1 ; should be replaced by the condition. 



r[{d/dt\oga{t)r + El^J-2<^l. 



(44) 
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Figure Captions 



Figure 1 < > obtained by numerical solution of Eq. (21) (randomly fluctu- 

ating off-diagonal) with Bq = 1, k = 0.05 for (a) A = 0, (h) 0.5, (c) A = 2.0 and 
(d) A = 10. The initial condition is ^i(O) = 1,^2(0) = 0. 

Figure 2 < |'0i(t)p > obtained by numerical solution of Eq. (24) (randomly fluctu- 
ating diagonal) with = 1, A; = 0.05 for (a) A = 0, (b) A = 0.5, (c) A = 2.0 and 
(d) A = 10. The initial condition is -01(0) = 1,-02(0) = 0. 

Figure 3 kAr as a function of Am^/ E for a randomly fluctuating magnetic field of 
rms value 1 x lO^C confined to a region of the solar convective zone, as described in 
the text. 

Figure 4 < {ipel"^ > for infinite distance as a function of Am'^/E. The lower(upper) 
curve around Am^/E = 1 x 10~^ is the curve for a noiseless (noisy) electron density. 
The rms value of the randomly fiuctuating noise is 2 % of the local density, sin^ 29 is 
changed from (a) 0.8, (b) 0.1, (c) 0.01 and (d) 0.001. The correlation length is given 
by Eq. (35) for each value of Am?/E. 

Figure 5 The same as Fig. (4) except that the rms value of the randomly fiuctuating 
density is 4 % of the local density. 
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